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ABSTRACT 

Aims. We search for the long-term variability of the y-ray sky in the energy range E > 1 GeV with 168 weeks of Fermi-LAT data. 
Methods. We perform a full sky blind search for regions with variable flux looking for deviations from uniformity. We bin the sky into 
12288 bins using Healpix package and use Kolmogorov-Smirnov test to compare weekly photon counts in each bin with a constant 
flux hypothesis. The weekly exposure of Fermi-LAT for each bin is calculated with the Fermi-LAT tools. We consider flux variations 
in the bin significant if statistical probability of uniformity is less than 4 x 10~ 6 , which corresponds to 0.05 false detections in the 
whole set. 

Results. We identified 117 variable sources, variability of 27 of which has not been reported before. Among the sources with 
previously unidentified variability there are 25 AGNs belonging to blazar class (11 BL Lacs and 14 FSRQs), one AGN of uncertain 
type and one pulsar PSR J0633+1746 (Geminga). The observed long term flux variability of Geminga has a statistical significance of 
S.lcr. 
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1. Introduction 

Time domain astronomy at different wavelengths from radio to 
very high-energy gamma-rays is developing very rapidly nowa- 
days. In the high energy (HE) range (> 100 MeV) great progress 
was a chieved with adven t of the gamma-ray telescope Fermi- 
LAT dAtwood et al.l 12009b . Its high sensitivity and almost uni- 
form sky coverage allow one to study variability of large number 
of sources at these energies on time scales from seconds to years. 
Sources that demonstrate the most variable behaviou r are active 
galacti c nuclei (AGN), primarily of the blazar type (Abd o et aD 
2010a). It is well known that these sources exhibit v ariability on 
\ differ e nt time scales and at different wavele ngths ( Carin i et al.l 
1 199 It lUlrich etal.1 Il997b IWelsh et al.l l201lt kani et all l2009t 
Bon ning et al.l l2009t ISoldi et all |2008[ iGhisellini & Tayecchiol 
20081: iRaiteri et al.ll2005b ICiprini et alTl2003t lUrrv et al.lll993l) . 
Studies of time variability of these sources are very important 
for better understanding of AGN engines; they are also essential 
for assessing the quality of spectral energy distributions obtained 
from multiwavelength observations made at different epochs. 

In this paper we perform a full sky blind search for variable 
sources. We bin the sky into equal area pixels and search for 
deviations of photon number counts from the uniformity in time. 



2. Data and method 

The LAT Pass 7 weekly all-sky data publicly available at Fermi 
mission websit^D were used in this work. The analysis covers 
the time period of 168 weeks from August 04, 2008 to October 
18, 2011, corresponding to mission elapsed time (MET) from 
239557417 s to 340622181 s. We use the 'Pass 7 Source' event 
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class photons with E > 1 GeV and impose an Earth relative 
zenith angle cut of 100° and rocking angle cut of 52°. 

We bin the dat a week by week using HEALPIX package 
dGorski et al.l l2005) into a map of resolution N s ide = 32 in galac- 
tic coordinates with 'RING pixel ordering. Total number of pix- 
els is equal to 12288 and the area of each pixel is 3.6 sq. deg, 
chosen according to the size of Fermi-LAT point-spread function 
(PSF) above 1 GeV which is approximately 1°. We estimate in- 
tegral weekly exposure for each pixel using the standard Fermi- 
LAT tools gtltcube and gtexpcube (ScienceTools-v9r23pl-fssc- 
20110726). 

For each pixel we count the number of photons in each of 
168 weeks and consider corresponding values of weekly ex- 
posure. Typically there are 2-10 photons in a pixel per week 
except bins with the brightest sources. Cumulative distribution 
functions (CDFs) Pit), £(?) for both photon counts and exposure 
are constructed. In the absence of variability, the photon counts 
would represent a random process with CDF proportional to £(f) 
and thus P(f) would follow £(f) with deviations caused by a fi- 
nite number of observed photons. Otherwise P(f) would not be 
statistically compatible with £(f). A Kolmogorov-Smirnov (KS) 
test is a natural and straightforward way to examine statistical 
compatibility of the observed photon counts with the distribution 
given by CDF £(f). The probability that both sets represent the 
same distribution could be estimated from the maximal value of 
distance between the functions fit) and £(f). An example CDFs 
for one of the pixels is shown in Figure[T]and the corresponding 
flux is shown in Figure [2] 

Implemented KS test is most sensitive to the variability at 
long scales (longer than a week), while transient bursts and flares 
at shorter time scales may be missed if they are not overwhelm- 
ingly strong. On the other hand, our method is sensitive to grad- 
ual moderate changes in photon fluxes without any prominent 
bursts, which could be missed by burst searching techniques. 
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Fig. 1. Plot of the cumulative functions G(t) and V(t) for pixel 
no.54 (/ = 261°, b = 82.°69). The difference could be easily 
seen and probability that the photon rate is constant is only P = 
4 x 1(T 80 . 
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Fig. 2. Flux of photons with energies larger than 1 GeV for pixel 
no.54 (4C+21.35, see Figure[T]i. The flux is in photons cirT 2 s _1 
units. 

KS probability is calculated for each pixel - we are interested 
in pixels with probabilities smaller than the threshold value Pq = 
4 x 10~ 6 . This threshold value is set to allow for penalty coming 
from the large number (/V p ; x = 12288) of trials: the detection 
criterion is chosen in such a way that the entire search would give 
a single false detection with the probability of Pq x N p \ x ~ 0.05. 
A map of the probabilities is presented in Figure|4] 

3. Results 

The total number of bins with probabilities smaller than the 
threshold value Pq is 151. The source identification for each of 
these bins is performed as follows. We consider the photons that 
arrived during several weeks at the epoch of the maximal flux. 
The center of mass of the spatial distribution of these photons is 
used as an initial estima te of the source location. We sear ch for 
sources from 2FGL (iThe Fermi-LAT Collaborationll201 lh cata- 
log in the circle with a radius equal to the error of center of 
mass estimation which is usually about 10' -15' for ~ 20 pho- 
tons). For 34 pixels no source is found with the initial esti- 
mate: in 3 1 of them the variability time pattern and photon spa- 
tial distribution leads us to identification of the source located 
in the neighbouring bin and already identified there. In three 
cases (pixels 2877 and 2749, PMN J1532-1315, and pixel 2299, 
SBS 0846+513) we have found no counterpart in 2FGL catalog 
and use Simbad astronomical database. Gamma-ray flares from 
the latter two sources (PMN J1532-1315 and SBS 0846+513) 
have been recently reported by Fermi-LAT ( Gas parrini & Cutinl 
1201 ltlPonato & Perkins 201 lh . but the activity period is outside 
of the time region of 2FGL catalog. In order to avoid false iden- 
tifications we performed an additional check in four cases when 



two different sources are residing in adjacent pixels: significant 
difference in the observed luminosity curves confirms that these 
detections are real. The results are presented in Table Q] 

The total number of ident ified sources is 117; variability 
of 55 of them was reporte d in dAbdo et al.ll2010at iTanaka etaD 
lIoTltlSchinzeT et al.ll20iTt) . 35 additional detections were made 
in numerous Astronomers 's telegrams (ATels, see Appendix) 
and on the Fermi-LAT blogQ. That leaves us with 27 sources for 
which the variability has not been previously detected (see Table 
[2] and Figures [5]-[9]l. We have explicitly checked that flux from 
these sources is not contaminated with the contribution from the 
Sun. 

For convenience we assign a variability type to the source: a 
gradual change in the photon flux is referred to as 'rate', if the 
whole variability is dominated by one or several flares we call 
it 'flare' and if these flares are observed by more or less pro- 
longed time span (typically, more than 20 weeks) we designate 
it as 'activity'. This morphological distinction is not totally un- 
ambiguous: several bright consequent flares could be defined as 
'activity'. On the other hand, gradual changes could take place 
on a time scales of several tens of weeks, thus fitting the 'activ- 
ity' type as well. 

The BL Lacs and the flat-spectrum radio quasars (FSRQs) 
are represented almost equally in the set of previously unidenti- 
fied sources: there are 1 1 BL Lacs, 14 FSRQs, and one AGN 
of uncertain type (PKS 0644-671). Also there is one pulsar 
(Geminga) in the list (see Table [2] and Figures [5]-[9]l. Both BL 
Lacs and FSRQs demonstrate two types of variability: gradual 
change of photon flux (10 out of 25) and flares or increased ac- 
tivity (15 out of 25). Variability of several sources was observed 
before in other energy ranges: a flare in the near infrared region 
was o bserved for B2 1732+38A (ATel #3504 dCarrasco et al.1 
1201 lh . 1 July 2011), VHE flare s of 1ES0806+ 524 were ob- 
served by MAGIC (ATel #3192 dMariottil l20l"Tb . 24 February 
2011), and EGRET obser ved a very bright flar e of PKS 2255- 
282 in December 1997 dMacomb et all Ill999h . Pulsed y-rays 
from the G eminga pulsar wer e observed with 1 year of Fermi- 
LAT data (AbdoetaDHoIoi), while the source were consid- 
ered non-variable. In this paper the long term flux change of 
Geminga pulsar is detected with KS probability of 2.3 x 10" 7 . 
In present study Geminga is the only pulsar demonstrating 
variability above our threshold. For comparison the KS prob- 
abilities for bins containing Vela, PSR J 1709-4429, and Crab 
pulsars are 25%, 78% and 1.2% respectively. In Crab case it 
was shown that the observed variability is cause d by processes 
in the Crab nebula rather than the pulsar itself (Buehle r et al.l 
1201 lh . Fermi-LAT collaboration also presented results of ob- 
servations of three other gamma-ray pulsars (J 1836+5925, PSR 
J202 1+4026, and PSR J0007+7303) where no flux variabil- 
ity was obs erved dAbdo et al.ll2010UlSaz Parkinson et alJl2Q10fc 
lAbdo et alJl20l2h . That makes the case of Geminga even more 
intriguing. 

We note that while 2FGL catalog contains 577 unidentified 
sources (out of 1873), 153 of which have flux higher than 2 X 
10~ 9 photons crrT 2 s _1 , none of them show variability above our 
threshold. 

It is also worth noting that sources not included in the Table|2] 
because of being reported either in ATels or on the Fermi blog 
could have the variability patterns that differ considerably from 
the reported one. As an example, the flare from the source MG2 
J 130304+2434 that took pl ace on 3 July 2009 ( week no. 56) 
was reported in ATel #2110 dHavs & Marellil2009h . On the other 
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Fig. 3. Variability of MG2 J 130304+2434. Change of rate could 
be easily seen. 

hand, Figure [3] shows that the flare occured during the high state 
of the source, with its flux slowly increasing from the start of 
the Fermi observations till approximately the 80-th week when 
it began to decrease. 

4. Conclusions 

A method for variable sources detection is proposed that uses 
the KS statistical test. The method is implemented for a full sky 
blind search for regions with variable flux at energies above 1 
GeV using Fermi-LAT 168 weeks data. The search leads to iden- 
tification of 1 17 variable sources, the variability of 27 of which 
has not been reported before. Among the sources with previously 
unidentified variability there are 25 AGNs belonging to blazar 
class (11 BL Lacs and 14 FSRQs), one AGN of uncertain type 
(PKS 0644-671), and one pulsar PSR J0633+1746 (Geminga). 
The observed long term flux variability of Geminga pulsar has a 
statistical significance of S.lcr. 

Appendix 

The following ATels are cited in t he text: ATel #1933 
dCorbel & Revesl [2009h . ATel #2048 dCiprinil l2009cb . ATel 
#2049 (ICiprinill2009al). ATel # 2104 (lLongo et al.ll2009h. ATe l 
#2110 riHavs & Mareliil l2009h ATel #2136 dCiprinil 120 09b). 
ATel #2243 dTanaka et al.ll2009l). ATel #23 16 dHavs & Escandd 
2009) . ATel #2402 dSokolovskv et ail l2010b. ATel #241 3 
dffill & Vandenbrouckd koiOb ATel #2539 dWallace] | 2010b . 
ATel #2583 dDonatol 12010). ATe l #2669 dCutinil I20T0B 
ATel #2783 dD'Ammandol l2010bh. A Tel #2829 dSchinzell 
120101). ATel #2860 dD'Ammandol l2010ah ATel #290 7 
dCannon & D'Ammandol koiOb. ATel #2943 dCiprinil I20T0I). 
ATel #3002 dP'Ammando & Vandenbrouckd2010b, ATel#3026 
dAllafort & D 7 A mmand ol2010b. A Tel #3171 dBuson & Basfierll 
120 1 lb. ATel # 3 1 92 dMariottil20 l ib. AT el #3207 dAUafortl201 lb. 
ATel #3452 dDonato & Perkinsll201 lh. ATel #34 45 dGasparrinil 
12011b. ATel #3504 dCarrasco et all 1201 lb. ATel #3579 
dGasparrini & Cutinil [201 lb. ATel # 3670 dSchinzel & Ciprrnil 
1201 lb . ATel #3793 dOiha et~aDl2Ql lb . 



work is performed at the cluster of the Theoretical Division of INR RAS. This 
research has made use of NASA's Astrophysics Data System and SIMBAD 
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Table 1. List of pixels demonstrating variability exceeding threshold value (P < 4 x lCT^). l,b are the galactic coordinates of the center of the 
pixel, N phot is the total number of photons observed in the pixel, 0_ 8 is the average flux from the pixel 0_ 8 = O/ICT 8 photons crrr 2 s -1 : the total 
number of photons divided by the total exposure, P is the KS probability, the designations of the identified source in the literature and in the 2FGL 
catalog are in the 8 th and 9th columns. Previo us refere nces to the variabilit y of the source are presented in the last column: P stands for paper 
dAbdo et al.ll201Qj) . P2 for dTanakaet alj201lh . P3 for dSchinzel et alll201ll) . ATNNNN for ATel #NNNN, and BNNN indicates that the outburst 
from the source was mentioned on the Fermi blog in the NNNth weekly report. ATNNNN with prefix VHE or IR indicate that flare was observed 
(and reported in corresponding ATel) in some other energy range: very high energy (larger than 100 GeV) or in the infrared. All the references for 
the ATels are listed in the Appendix. 
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Table 2. List of sources with previously unreported variability. Additional columns describe source type: BL Lac (B), FSRQ (Q), 
AGN of uncertain type (AGN), or pulsar (PSR), variability type: gradual increase or decrease in photon flux rate (R), flares (F) or 
longer period of increased activity (A); in case of flares or activity the temporal localization of events is given in the last column. 



* 



-40.0 



0.00 



Fig. 4. Map of Fermi-LAT variability at 1 GeV in galactic coordinates. Pixel color represents base 10 logarithm of Kolmogorov- 
Smirnov probability of the uniformity of the observed flux. The galactic center is in the center of the figure, I = 180° is on the 
left. 
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Fig. 5. Luminosity curves for variable sources listed in the Table[2] The flux is in photons cm 2 s 1 units. 
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Fig. 6. The same as in Figure[51 
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Fig. 7. The same as in Figure[5] 
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Fig. 8. The same as in FigureO 
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Fig. 9. The same as in Figure|5] 
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